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An unidentified 3.5 keV line from X-ray observations of galaxy clusters lias been reported re¬ 
cently. Although still under scrutiny, decaying dark matter could be responsible for this signal. We 
investigate whether an axino with a mass of 7 keV could explain the line, keeping the discussion 
as model independent as possible. We point out several obstacles, which were overlooked in the 
literature, and which make the axino an unlikely candidate. The only viable scenario predicts a 
light metastable neutralino, with a mass between 0.1 and 10 GeV and a lifetime between 10~ 3 and 
10 4 s. 


I. INTRODUCTION 

There has been interest recently in an unidentified 3.5 
keV line in X-ray observations of galaxy clusters mm- 
Despite its physical origin being subject to debate El¬ 
ls], there is still some room to speculate that decaying 
dark matter is responsible for this signal. The first ob¬ 
vious dark matter (DM) candidate in this context is a 7 
keV sterile neutrino [I] , but many other alternatives have 
been proposed. Some authors have pointed out that a de¬ 
caying axino could explain the line EMI- In this paper 
we examine carefully the conditions under which a 7 keV 
axino would produce the observed X-ray line. Due to the 
large number of parameters at disposal in supersymmet¬ 
ric models, it is hard to exclude with certainty the axino 
scenario. However we show that various constraints leave 
almost no room available in the parameter space of these 
models. Therefore we deem the axino an unlikely can¬ 
didate to explain the line. The authors of Ref. m also 
mentioned the axino as an unnatural explanation of the 
X-ray line. In this work we elaborate on the physical 
arguments that lead to such a conclusion. 

There are a few reasons why it is appealing to con¬ 
sider a model with a light axino. First, introducing the 
axion multiplet (axion, axino, saxion) in models of su¬ 
persymmetry (SUSY) solves the strong CP problem [13] . 
Second, the axion and the axino can both be DM candi¬ 
dates. Third, if the axino has a mass in the keV range, 
it is warm DM and it could help reconcile some small- 
scale structure issues HMEZ] of cold DM. In the context 
of R-parity violating (RPV) SUSY [T8|, a light axino is 
unstable. Its lifetime is still longer than the age of the 
universe, but it can decay into a neutrino and a photon. 
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It is this channel that would produce the 3.5 keV line. 
RPV models also have the virtue of explaining the null- 
results for SUSY searches at the LHC without introduc¬ 
ing a little hierarchy problem: the limits on the masses 
of the sparticles become much weaker [191123] . The axino 
could also produce an X-ray photon in R-parity conserv¬ 
ing SUSY. This is possible if the lightest neutralino is 
very light or even massless. The axino can then decay 
into the neutralino plus a photon [12] . We comment also 
on this possibility, below. 

This letter is organized as follows. In sec. [TT] we review 
how the axino abundance depends on the reheating tem¬ 
perature, in sec. m we discuss various constraints that 
the 3.5 keV line puts on models with a decaying axino. 
We conclude in sec. [IV] 


II. AXINO RELIC DENSITY 

Any supersymmetric model with an axino contains also 
an axion. The latter, as an invisible axion, is always a 
good DM candidate, while the former is suitable only if 
it is sufficiently long lived. We are here interested in a 
scenario where the axino constitutes almost the entire 
DM budget of the universe. This is typically realized for 
low values of the axion decay constant ( f a ~ lO 10 GeV), 
in which case the axion DM component is suppressed and 
can be neglected. See e.g. Ref. [M] for a review on axion 
DM. 

The axino as a DM candidate has been widely explored 
in the literature ]25ti32j . There are two axion/axino sce¬ 
narios of interest: the Dine-Fischler-Srednicki-Zhitnitsky 
(DFSZ) model [33) [33], and the Kim-Shifman-Vainshtein- 
Zakharov (KSVZ) model |35l [35] . As far as the axino 
is concerned, the important difference between the two 
lies in the mass, M$, of the heaviest Peccei-Quinn (PQ) 
charged and gauge-charged matter supermultiplet in the 
model [31] • In the DFSZ case is the Higgsino mass, of 
order < TeV, while in the KSVZ case it is the mass of the 
heavy vectorlike quarks, typically of order f a > 10 9 GeV. 
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This difference leads to two different relic axino abun¬ 
dances as well as to a different dependence on the reheat¬ 
ing temperature, 7 rh- For DFSZ models we have |29j 



which does not depend on Trh, as long as Trh is larger 
than 1 roughly 1 TeV [52] . The relic density drops very 
quickly for lower values of Trh (see e.g. Fig. 3 in 
Ref. j2£i]). From eq. ([l]) we see that for a 7 keV DFSZ ax¬ 
ino with f a < 10 10 GeV the reheating temperature has to 
be below 1 TeV in order to avoid overabundance. Recall, 
obervationally Hdm^ 2 = 0.1199 ± 0.0027 137]. 

For the KSVZ case 29] : 



which is valid for 1 TeV < Trh < ~ f a ■ Also in 

this case the abundance drops very quickly for reheating 
temperatures below about 1 TeV. 


A. Including a light gravitino 


Without tuning and/or cancellations, the axino is ex¬ 
pected to be heavier than the gravitino [38] [39]. This 
puts further constraints on our scenario. First, the axino 
can decay into a gravitino and an axion with a lifetime 

m 


>G+a 


~ 6 - 10 


28 



(3) 


Here m 3 / 2 denotes the gravitino mass. If the 7 keV ax¬ 
ino is to explain the X-ray line we must require its life¬ 
time to be roughly greater than the age of the universe, 
T a^G+a > 10 18 s > otherwise it would have decayed away. 
This implies m 3 / 2 > 10 _5 keV. A gravitino lighter than 
roughly 100 eV is a hot dark matter candidate. Its relic 
abundance 00 ], 

fig DM /i 2 - 0.1 , for m 3/2 < 100 eV, (4) 

is constrained to be less than 3% of the total DM abun¬ 
dance [5T], which implies m 3 / 2 < leV. 

In the range 1 keV < m 3 / 2 < 7 keV, the gravitino is a 
warm dark matter candidate and its abundance depends 


1 The exact value of Trh depends on the SUSY spectrum, but it 

is expected to lie in the TeV range. 


on Trh 02]: 

^ = 0 ' 27 (loStSv) O • for "*^ > lkcV - 

(5) 

In this case we require a reheating temperature below 
1 — 10 GeV so that the gravitino contribution to the DM 
density is much smaller than that of the axino. However 
for such low values of Trh the axino relic density, even 
with f a ~ 10 9 GeV, is highly suppressed, fl^h 2 <C 0.12. 
Therefore, we exclude this case. 

Thus if the axino constitutes most of the DM and pro¬ 
duces the 3.5 keV X-ray line, the gravitino mass is re¬ 
stricted to the small window 

10” 2 eV < m 3/2 < 1 eV. (6) 

As both the axino and the gravitino are in the keV 
range and below, one might worry about their contribu¬ 
tion to the relativistic degrees of freedom at Big Bang 
Nucleosynthesis (BBN). This turns out not to be a prob¬ 
lem in most cases. When the universe reheats to 10 - 
100 GeV, both the axino and the gravitino are out of 
thermal equilibrium [29] . The subsequent annihilations 
of SM particles heat up the photon bath so that the pho¬ 
ton temperature at BBN is higher than the respective 
temperatures of the axino and the gravitino. As a conse¬ 
quence their contribution to N e g is well within the bound. 
If there is also a very light or massless neutralino in the 
spectrum one has to worry about the constraint from 
N e g , as we discuss later. 


III. THE AXINO AND THE X-RAY LINE 

To explain the 3.5 keV X-ray line via the decay of a 
DM particle X, one needs a decay rate of Fa-->. 7 +... ~ 
(10 28 s)" 1 ~ 10 -53 GeV [T], assuming that the decay¬ 
ing DM constitutes 100% to the relic abundance. If the 
decaying DM is a fraction k < 1 of the total DM, then 
the corresponding decay rate has to increase by l/k to 
explain the signal. 

There are two scenarios for a 7 keV axino to produce 
the 3.5 keV X-ray line, where R-parity is respectively 
conserved or violated. In either case the starting point is 
the following Lagrangian for the coupling of the axino to 
the SM gauge bosons and their gaugino supersymmetric 
partners: 

£aAV = i(,J 2 y CaBBBB fJ,l ' + 

9 2 2 C aWW W a W a ^ + 5 2 G“G“^) . (7) 

Here a is the axino mass eigenstate, while the gauginos 
and gauge fields are gauge eigenstates; g r are the gauge 
couplings. These interactions are the supersymmetric 
version of those of the axion with gauge fields in the 
low energy (non-supersymmetric) Lagrangian [43]. The 
coefficients G are model dependent and are typically of 
order one [22] , 
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A. R-parity conserving SUSY 

In the R-parity conserved case, if the bino is lighter 
than the axino, the X-ray line could be produced by the 
decay a —> B + 7 . This was pointed out in Ref. [T2] ■ The 
axino partial decay rate with a massless bino is 



(9) 


Note, a massless bino is consistent with all data provided 
the sfermions are heavy enough [151 - HT] . To match the 
rate needed for the X-ray line this scenario requires f a > 
10 14 GeV. This immediately excludes the DFSZ axino, 
whose relic abundance would be much too low. 

From eq. ([ 2 ]), a KSVZ axino with Trh ~ 10 12 GeV 
would seem viable. However this scenario is strongly 
disfavored by two arguments. First, the abundance of 
axions produced via the misalignment mechanism would 
be far too high [24] for f a ~ 10 14 GeV, unless one tunes 
the initial misalignment angle to very small values. Sec¬ 
ond, a massless bino together with a light gravitino would 
contribute to the relativistic degrees of freedom [45] , giv¬ 
ing a value A N e g > 1, in strong tension with BBN and 
with the data from the Cosmic Microwave Background 
(CMB) g5]. 


B. R-parity violating SUSY 


Here r v .g (r s) parametrizes the mixing of the neu¬ 
trino mass eigenstate with the gaugino gauge eigenstate 
B (W°). The lifetime of the axino to explain the X-ray 
line requires r 2 ^ (r 2 ^) t° be of order 10~ 12 for f a fixed 

at its lowest possible value [51], 10 9 GeV. 

One the outstanding features of RPV models with lep¬ 
ton number violation, is that they automatically pro¬ 
vide for massive neutrinos. Assuming that the neutrino 
masses solely arise from the RPV sector, we can estimate 
bounds on the mixings G/(b vv) as follows. Neglecting 
loop contributions, the terms eiLiH u lead to one massive 
neutrino |521l54l and two non-vanishing lepton mixing an¬ 
gles, which we take to be 9 13 and d 2 3 [52] • The neutrino 
mass is given in terms of the model parameters as 


ffjMi + g\M 2 . y , 2 
4det A4 1 1 ’ 


Xo 


( 12 ) 


where [52| 


A i = nidi — Vd,£i , i = e,n,r (13) 

are the alignment parameters and tc, = (p.j) are the sneu- 
trino VEVs. Furthermore det A4 Xo denotes the determi¬ 
nant of the 4x4 neutralino sub-mass-matrix of the MSSM 

det M Xo = -n 2 M 1 M 2 + ~nv u v d (g 2 M 1 + Si M 2 ) , (14) 

where v u = (H u ), Vd = (H d ) are the Higgs VEVs. The 
parameters A, are related to the two remaining neutrino 
mixing angles [52]: 


In the context of R-parity violation, the terms 
eiLiH u , i = e,n,T hr the superpotential introduce mix¬ 
ing among the neutrinos and the Higgsinos. The modified 
scalar potential also results in non-zero sneutrino vac¬ 
uum expectation values (VEVs), which introduce mix¬ 
ing between the neutrinos and the bino and the neutral 
wino, respectively. In this case the RPC axino decay 
in sec. |III A| automatically includes the decay channel 
a —> Pi+ 7 . The new partial decay rate is simply modified 
by the appropriate mixing angles 2 



2 Note this decay would also occur with pure trilinear R-parity 

violation via the resulting sneutrino vev’s El 


tm(,1 3 = - (A;+%),/• ■ = (15) 

These angles are measured [55j to be ~ 7r/20 and 
$23 ~ 7t/ 4, thus we have A M = A r = A and A e = 0.23 A. 
The cosmological bound on the sum of the neutrino 
masses [49], JV m Vi < 0.23 eV, in our case amounts to a 
bound on the single massive neutrino, and thus 

A < (3.2 x 10- 13 TeV) 1 / 2 ( ^■ (16) 

v ; \g 2 2 M 1+ glM 2 J y J 

From Refs. [52] [SJ we can work out the mixings analyt¬ 
ically, in terms of the supersymmetric parameters. For 
instance 

_ _ giM 2 A M A r /r _ giM 2 Ag nr> 

Tv2B ^A2+A 2 - det M X0 ~ V2det M Xo ’ 


and similar expressions for the other mixings, 
bound on A translates into 

Then the 

r 2 ^ < (1.1 x 10" 14 TeV)A4 _1 

(18) 

r^~ < (4.6 x 10- 13 TeV)A4 _1 

(19) 

rl 3§ < (2.8 x 10" 16 TeV)A4 _1 , 

(20) 
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with 


M~ x 


__ 

(ffl-Mi + 9iM 2 ) det M Xo ’ 


( 21 ) 


and similar bounds for r 2 —. Given the null SUSY 
searches at the LHC so far, it is reasonable to expect 
the parameters Mi,M 2 ,ijl to be of order TeV or larger, 
and det A4 Xo ~ TeV 4 , in absence of cancellations. In this 
case the bounds simplify to 


.2 

vi B 

< 

2.3 x 10" 15 

.2 

vi B 

< 

9.2 x 10" 14 

.2 

v 3 B 

< 

5.6 x 10~ 17 


( 22 ) 

(23) 

(24) 


These mixings are very small and thus the resulting axino 
decay rate too slow to explain the X-ray line. 

Perhaps the assumption that all the SUSY parameters 
and masses are at least around a TeV is too strict. Sup¬ 
pose that the lightest neutralino, Xii has a mass of order 
GeV and the other neutralinos are at the TeV scale. Then 
detA4 Vn ~ 10 -3 TeV 4 and the bound on r 2 - becomes 

AU 1/2-D 

of order 10 -10 , which is enough to fit the line. However 
this scenario faces another problem. A GeV neutralino 
has two decay channels: one into an axino and a photon 
with 




1 


1287T 3 /' 


1 

aBB 


fj:) cos2 ° w ’ 


>2+7 


60 s 


GeV 


fa 


10 9 GeV 


(25) 


the other into a neutrino plus leptons (or pions) via an 
off-shell Z boson, with [56 


- x?-vW+2- 


<b q2 m xg 

1024tt 3 Mf ’ 



GeV 


m v o 

Xi 


(26) 


For m x o < 1 GeV the dominant decay mode is into ax¬ 
ino plus photon, while for ?n x o > 1 GeV it is into neu¬ 
trino plus leptons. These decay rates are low enough 
that the neutralino freezes out before decaying. As a 
weakly interacting particle (WIMP), at freeze-out it has 
roughly the right DM relic abundance, fl x o h 2 ~ 0.1. If 
Xi is mostly bino, it will be slightly overabundant, while 
if it is mostly wino or higgsino it will be underabun- 
dant m • In either case our subsequent conclusions do 
not change. When y 9 decays it produces energetic pho¬ 
tons, which are subject to constraints from BBN and 
CMB. The window 10 MeV < m v o < 100 MeV, corre- 
sponding to 3 x 10 4 s < r x o < 3 x 10 7 s, is excluded by 
bounds from photodestruction of D and photoproduction 
of D+ 3 He 21] ■ The window 300 keV < m x o < 10 MeV, 


corresponding to 3 x 10' s < r x j < 10 12 s, is excluded 
by bounds from CMB spectrum distortions [SHI S51 . 
The window 7 keV < o < 300 keV, with a lifetime 

A1 

r x o > 10 12 s, is excluded by CMB constraints on late de¬ 
caying particles [60] . For nn x o < 7 keV the neutralino is 
lighter than the axino and we are back to the situation of 
the R-parity conserving scenario, which is strongly disfa¬ 
vored as we discussed in the previous section. 

The window 100 MeV < m x o < 10 GeV, correspond¬ 
ing to 10~ 3 s < r x o < 3 x 10 4 s, cannot be as easily ex¬ 
cluded, so it is viatle in principle. In this range most of 
the neutralinos decay around BBN time, which could be 
in some tension with the success of BBN itself. However 
a detailed study, beyond the scope of the current work, 
is needed to come to a definite conclusion. For m v o >10 
GeV one quickly hits the bound from neutrino masses. 

We conclude that also the RPV scenario is strongly 
disfavored. 


IV. SUMMARY 

Motivated by recent claims of detection of an X-ray 
line at 3.5 keV we have investigated whether a decay¬ 
ing axino with a 7 keV mass could explain the signal. 
The R-parity conserving scenario is strongly disfavored 
mostly because it requires a very light bino, which to¬ 
gether with a light gravitino would contribute to the rel¬ 
ativistic degrees of freedom, AN e s > 1, in contradiction 
with BBN and CMB bounds. RPV models require a 
mixing neutrino-bino or neutrino-wino which is typically 
too large and excluded by the cosmological bound on 
neutrino masses. To evade this bound one is forced to 
take the lightest neutralino with a mass around a GeV 
or below. Such a neutralino is long-lived and decays into 
energetic photons. If it is lighter than 100 MeV the sce¬ 
nario is excluded by BBN and CMB constraints. These 
arguments are generic, model independent, and make the 
7 keV axino a very unlikely candidate for the observed 
3.5 keV line. 

If the lightest neutralino has a mass between 100 MeV 
and 10 GeV, the RPV scenario is still viable. The cor¬ 
responding lifetime of the neutralino ranges from r x o = 
0(lO _3 s), for m x o = 10 GeV, to r x o = C(10 4 s), for 
m x o = 100 MeV. This is effectively stable for collider 
physics. Such a neutralino is very hard to observe in the 
laboratory, with no present bounds [46] . The low mass 
range is excluded by supernova cooling if the selectron 
is lighter than about 500 GeV g5j. The proposed SHiP 
facility at CERN is most likely not sensitive to these life¬ 
times, due to the restricted geometry [6L. A direct mea¬ 
surement at an e + e _ linear collider is also unlikely [62] , 
This range of lifetimes is too short for possible astrophys- 
ical signatures. Therefore it is hard to verify whether this 
particular scenario is realized or excluded. 

It might still be possible, at the expense of fine tun¬ 
ing, to find other corners of parameter space where the 
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axino has the right abundance and decay rate to fit the 
line. One would have to do a numerical scan for a spe¬ 
cific model, which is beyond the scope of this work. If 
such a corner were found, our arguments suggest that it 
would have a low axion decay constant, f a ~ 10 9 GeV, at 
the edge of the astrophysical bound, and a low reheating 
temperature, below a TeV (for a KSVZ model) or even 
around a few tens of GeV (for a DFSZ model). This is 
troublesome for most of the proposed baryogenesis mech¬ 
anisms. 


Note added 

As this work was being completed, Ref. [63] appeared, 
which casts further doubts on the DM interpretation of 


the 3.5 keV line. 
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